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ABSTRACT 



Measurements of forward scattered light were used to 
determine the particulate behavior in the combustor and 
across the exhaust nozzle of a small solid propellant rocket 
motor. A MALVERN 2600 was successful in measuring 

multimodal particle distributions, and a locally designed 
system was successfully applied to monomodal distributions. 
The propellant tested was a reduced smoke propellant 
containing one-percent zirconium carbide by weight. The 
results suggested that ZrC reacts to completion on or very 
near the propellant surface. Zirconium oxide (Zr0 2 ) 
agglomerates are believed to form on the surface. Surface 
agglomerate size decreased with increasing pressures to 
approximately 500 psig. Agglomerates broke up in the nozzle 
convergence to approximately 18 microns, independent of the 
pressure or nozzle entrance particle size.- Particle 
breakup occurred in the exhaust plume.. System limitations, 
sources of error and suggested improvements are discussed. 
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I. 



INTRODUCTION 



The use of metallized propellants offers the advantage 
of improved solid rocket motor performance in terms of 
increased energy, specific impulse, and suppression of high 
frequency acoustic combustion instabilities. Aluminum is 
currently the primary metallic additive of choice when 
exhaust plume visibility is not a mission consideration. 
However, the requirement for reduced smoke propellants for 
use in tactical air-to-air applications has dictated the use 
of alternate high energy propellants which utilize only 
small amounts of additives such as zirconium, zirconium 
carbide, and aluminum oxide to provide acoustic 
stabilization . 

The improvements achieved via such additives are not 
without penalty, certain loss mechanisms are introduced. 
The specific impulse efficiency of a propellant containing 
metal is generally less than that of a base propellant due 
to incomplete combustion of the metal in the motor and the 
presence of condensed metal and metal oxides in the nozzle 
flow. Particle-gas velocity lag represents the primary 
source of performance loss within the exhaust nozzle. 
Thermal lag generally has less effect on performance. Yet, 
neglecting thermal and velocity lag as loss mechanisms, the 
simple fact that the condensed particulates do not expand as 
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the other combustion gases expand, results in a reduction of 
available specific impulse. 

Netzer et al. [Ref. 1], point out that the behavior of 
metal as it leaves the burning surface is not a well-defined 
phenomenon. Some particles upon reaching the burning 
surface depart immediately, while others agglomerate on the 
surface before passing into the gas flow. The mechanics of 
metal droplet combustion as discussed by Schorr [Ref. 2], 
indicate that resultant particles in the range of 1-2 
microns and smaller are of significance in determining 
exhaust plume signature, while those particulates of 5 
microns and greater are important in the determination of 
two-phase flow losses. The problem of predicting 
performance depends on an accurate determination of these 
two-phase flow losses which are a strong function of 
particle size. The aluminum oxide that results from the gas 
phase combustion of aluminum is generally less than two 
microns in size. These small particles have little velocity 
lag. Thus, making them smaller would serve little benefit. 

Several complex computer codes attempt to model the 
combustion process in order to predict solid rocket motor 
performance. The SPP (Solid Performance Program) [Refs. 
3,4], is one such code. These programs are semi-empirical 
in nature and require an a priori knowledge of particle size 
distribution entering the exhaust nozzle. This distribution 
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must be experimentally determined, as the problem is 
currently intractable with regard to analytic methods. 

The following experimental methods have been employed in 
previous work at the Naval Postgraduate School: 

(a) collection of exhaust products 

(b) light transmission measurements 

(c) light scattering measurements 

(d) holography and high speed motion pictures. 

A summary of previous research conducted at the Naval 
Postgraduate School combustion laboratory is contained in 
Reference 1. The present investigation focused on the use 
of light scattering techniques to determine particle size. 
The method is non-intrusive and suited to measurements in 
both the motor and exhaust nozzle environments. The speed 
of the flow precludes successful application of most 
holographic techniques within the exhaust nozzle. 

Heretofore most studies have been directed at aluminized 
composites, and much is known about the combustion of 
aluminum. The behavior of other metallic additives have not 
been as fully explored. For example, zirconium carbide is 
used as a stability additive. However, it is not known 

whether surface agglomeration occurs, or whether the 

zirconium burns to zirconium oxide. Both of these processes 
will affect the size distribution and density of the 

particles in the gas phase. These in turn affect the 
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ability of the particles to provide damping for a specific 
frequency of pressure oscillation. 

The goals of this investigation were: 

(1) to refine previously developed experimental and 
diagnostic techniques, and 

(2) to use the apparatus to determine the effects 
operating pressure and position within the motor and 
exhaust nozzle on the particle size distributions 
from a reduced smoke propellant. 

The propellant studied was a low smoke, composite which 
contained one percent zirconium carbide by weight. The 
behavior of particulates at various positions in the motor 
and exhaust nozzle was examined as a function of chamber 
pressure . 
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II . THEORETICAL BACKGROUND 



The performance of a solid rocket motor utilizing 
metallized propellants is significantly dependent on the 
diameter of the condensed particulates present in the flow. 
A non-intrusive , in-situ method is the most desirable means 
to determine particle size, since the measurement is real 
time and the environmental conditions (pressure, 
temperature, gas composition, etc.) are undisturbed. The 
method chosen for the present research was the dif f ractively 
scattered light technique. Very large particles cannot be 
measured with this technique due to the transmitted light 
interference with the light which is scattered at very small 
forward angles. Very small particles scatter at larger 
angles, but the intensity is low, making detection 
difficult. Absorbtion and refraction also become important 
for very small particles. Particle concentrations are also 
limited to those which do not exhibit significant multiple 
scattering effects. Within these limits, the mean particle 
size is readily determined from diffraction techniques, but 
size distributions cannot easily be determined. In general, 
a distribution type must be selected a priori, i.e., Rosin- 
Rammler, log-normal or other such distribution must be 
assumed. Yet, the diffraction technique is non-intrusive 
and can be used in the internal motor and nozzle 
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environment, where velocities and particle concentrations 
can be high. 

The scattering and absorbtion of a plane monochromatic 
electromagnetic wave by an isotropic homogeneous sphere is 
described by Mie's solution to Maxwell's equations [Ref. 5], 
and the general theory of scattering developed by Mie is 
presented by Van de Hulst [Ref. 6]. According to Mie's 
theory, the characteristics of scattered light depend upon 
the particle size, shape, internal structure, refractive 
index with respect to the surrounding medium, the scattering 
angle and the incident wave polarization and wavelength 
[Ref. 7], The solution to Maxwell's equations for multiple 
scattering in a concentrated sample, where interaction of 
scattering centers is significant, presently represents an 
intractable problem, even if boundary conditions could be 
specified. Even the solution without multiple scattering is 
quite complex and difficult to incorporate into diagnostic 
equipment. 

When the particles are much larger or much smaller than 
the wavelength of the incident light, the equations simplify 
to Fraunhofer diffraction and Rayleigh scattering, 
respectively. The light scattered from particles within 
these restricted ranges can be measured and used to readily 
determine particle size. The size of particles of interest 
in solid propellant combustion depend on the application. 
When considering two phase flow losses, particle sizes of 
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interest are generally five microns or greater. Jones [Ref. 
8], has shown that Fraunhofer diffraction can be used for 
particles as small as one micron when a He-Ne laser is used. 
Fraunhofer diffraction is applicable in the case of larger 
particles because the forward scattered light in the center 
lobe is due principally to diffraction. The mechanisms of 
reflection and absorbtion do not significantly contribute to 
the center lobe in this instance (this was shown by 
Hodkinson [Ref. 9]). If these latter effects are present to 
a significant extent, then the complete Mie equations must 
be used. 

The classical results of Fraunhofer diffraction from 
particles, i.e., the discrete rings associated with a 
monodispersion, are typically not observed in practical 
combustion experiments. Typically, the particulates are 
polydispersed and, in many cases, a multimodal 
polydispersion is present. The intensity distribution due 
to a monodispersion is given by Fraunhofer diffraction [Ref. 
10] as, 



1 ( 6 ) 

1 ( 0 ) 



2J (a sin 0) 2 

r ± ] 

1 a sin 0 



(1) 



where is the first order Bessel function, a = ^ D/X , 9 is 
the scattering angle for the forward direction, D is the 
particle diameter, and \ is the wavelength of the incident 
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light. When considering small angles (0 < 15°), the 

expression simplifies to, 



1 ( 0 ) 

KO) 



2 J. (a0) 2 



(2) 



Dobbins [Ref. 11] showed that the normalized integrated 
intensity of forward scattered light from a polydispersion 
of large particles is given as. 



1 ( 0 ) 

KO) 



D °° 2J (a0) 2 

J [ -(^0)— 1 1 N r (D)DdD 



/ N (D)D 4 dD 
0 r 



( 3 ) 



where N r (D) is a distribution function describing the 
proportion of particles with diameter D in the sample. 
Dobbins et al. [Ref. 11] used the Upper Limit Distribution 
Function (ULDF) developed by Mugele and Evans [Ref. 12] and 
the same approach was followed in a portion of this 
investigation. The volume-to-surface mean diameter (Sauter 
mean diameter) D 32 of a polydispersion is defined by, 

D 

00 3 

/ N (D)D dD 

n - 0 r 

°32 = — d * 

00 

/ N (D)D 2 dD 
0 r 
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The significance 


of D 32 


(versus D 4 - 3 , etc.) is 


that the 


illumination profiles for 


different 


size distributions of 


the ULDF plotted 


against 


OtD 32 0 /A) 


are brought 


into near 


coincidence [Ref. 


11 ] . 








Buchele [Ref. 


10 ] showed that 


for tt D 32 G/ X 


less than 



three, a gaussian curve can be used for the scattered light 
intensity which closely matches the theoretical intensity 
profile for a polydispersion obtained by integrating 
Equation (3). The gaussian expression presented by Buchele 
is given as, 



= exp[-(.57^D 32 G/X)2 ] 



(5) 



The data reduction technique developed by Netzer and Harris 
[Refs. 1,13] is based on Buchele' s work. The method 
involves comparison of the experimentally obtained intensity 
profile with that produced by Equation (5) . Use of Equation 
(5) requires the value of 1(0), the forward scattered 
intensity at theta equal zero. Since transmitted light 
dominates at very low scattering angles, 1 ( 0 ) cannot be 
measured. A ratioing technique eliminates this problem. 
Equation (5) can be used to ratio I(0 2 )/I(0i) as follows, 



I (6 2 ) 
1 ( 9 ].) 



exp{-D 32 ( 



121,2 

\ > 



( 0 2 2 




( 6 ) 
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This eliminates the need to determine 1(0). However, I(5]_) 
at the "minimum" forward angle ^ must be taken from the 
experimentally measured profile at a location where the 
transmitted light effects are not significantly present. 
The resultant profile I(9 2 )/I(Qi) versus 0 2 can then be 
generated for a specified D 32 . D 32 is varied interactively 

until the profile generated by Equation (6) matches the 
measured profile. This procedure yields D 32 . However, this 
approach has an obvious disadvantage in that it is only 
accurate if the polydispersion is monomodal. 

A Malvern 2600 particle sizer was also used in the 
present investigation. The Malvern software is similarly 
based upon Fraunhofer diffraction. However, it is based on 
the inversion of the measured scattered data to obtain a 
particle size distribution. The Malvern allows use of a 
selected type of particle size distribution (Rosin Rammler, 
log-normal, normal) or of a model independent mode. The 
model independent method is described by the Malvern 2600 
Reference Manual [Ref. 14] as a constrained non-linear least 
squares method, which is able to measure multi-modal 
particle distributions with high resolution. The obvious 
advantages of the Malvern 2600 are the ability to measure 
multi-modal particle size distributions, and it has been 
independently tested. Dodge et al. [Ref. 15] conducted a 
comparison of the Malvern device and a phase doppler system 
and found that there was excellent agreement in point 
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measured Sauter mean diameters between the deconvoluted 
laser diffraction data and the phase doppler data. The 
Malvern system is currently one of the most widely used 
laser diffraction devices and it shows great promise for 
application to particle sizing in solid rocket motors. 

Other work has been done in the area of analytic methods 
that appears to provide the basis for refinements in 
particle sizing and determination of particle size 
distributions. Koo and Hirleman [Ref. 16] evaluated the use 
of integral transform techniques on Fraunhofer diffraction 
patterns. These techniques offer the advantage of direct 
inversion of scattering data to obtain a size distribution. 
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III. EXPERIMENTAL APPARATUS 



A. INTRODUCTION 

Laser diffraction systems have become widely used in the 
field of particle and droplet size measurements. Fraunhofer 
diffraction is the underlying principle involved. Hodkinson 
[Ref. 9] demonstrated that for particles greater than 
several microns in size, the primary contribution in the 
near forward small angle region is from diffraction rather 
than reflection or refraction and Jones [Ref. 8] has proved 
that in this same region, Fraunhofer diffraction may be 
applied to particles as small as one micron when a helium- 
neon laser is chosen as the source of illumination. Figure 
3.1 details a simple schematic for a laser light diffraction 
particle sizer. Such a system typically consists of a He-Ne 
laser source, beam expander (spatial filter and collimating 
lens), Fourier transform lens, and a detector. 

Typically, small particles diffract light at large 
angles and large particles diffract light at small angles. 
Figure 3.2 depicts the diffraction process. Note that the 
undiffracted light is focused to a single point. 

B . ROCKET MOTOR 

The series of experiments conducted required the use of 
two small solid propellant rocket motors. The motor 
depicted in Figure 3.3 is the same as that used by Pruitt 
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Figure 3.1 Schematic of a Laser Diffraction System [Ref. 16 ] 



Detector measures 
Integral scattering of all 




Figure 3.2 



Scattering Properties of Particles as 
a Function of Size [Ref. 14] 



14 




Figure 3.3 Motor Components 




Figure 3.4 Short Motor Components 
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[Ref. 17]. The chamber diameter was 2.00 inches and the 
propellant was cylindrical, 1.98 inches in diameter and 
1 . 00-1 . 25. -inches thick. All tests conducted were done with 
end burning grains. The long motor (Fig. 3.3) was used to 
obtain in-situ measurements of particle size prior to the 
nozzle entrance and in the exhaust. Figure 3.4 depicts the 
short motor which was used to obtain in-situ particle size 
measurements just aft of the propellant burning surface. 
This particular motor had been used previously in connection 
with holographic imaging. For additional specifications see 
Reference 18. 

Both motors were windowed, and the windows were kept 
clean using a nitrogen purge. The nitrogen was "spread" 
using a sintered metal sleeve which was inserted into the 
aperture. 

Ignition was accomplished using a BKN0 3 ignitor which 
was fired by means of a nichrome filament energized by a 12 
VDC source. 

The propellant was bonded into the motor casing with a 
self-vulcanizing silicone rubber (RTV) . The propellant 
studied was a reduced smoke propellant with 87% solids 
loading. Table 3.1 presents the major propellant ingredient 
weight percentages. 

C. LASER LIGHT SCATTERING APPARATUS 

The locally designed laser diffraction system makes use 
of an 8 mWatt He-Ne source in the same general configuration 
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TABLE 3 . 1 



PROPELLANT INGREDIENTS 



Ingredient 



Approx. Wt. % 



R4 5M 



9.2 



RDX 



4.0 



AP 



82.0 



ZrC 



1.0 



depicted in Figure 3.1. (A complete description of the 
apparatus is contained in Reference 20 [see Figures. 3.5 and 
3.6]). The scattered light is passed through a wide-pass 
laser line filter and a 500 mm focal length Fourier trans- 
form lens, which focuses the light onto the detector. The 
detector is a linear photodiode array (Reticon RL 1024G) 
which has 1024 photodiodes on 25 micron centers. The array 
is mounted on a "mother" board which contains the sample- 
and-hold circuits. Each photodiode measures the intensity 
of the scattered light and outputs a proportional voltage. 
The array is scanned at a rate of 33 kHz and a 6 msec blank- 
ing pulse is used to sequence the sweeps. The range of 
scattering angles within the field of view is determined by 
a number of factors. The upper limit of scattering angle is 
determined by the diameter and focal length of the focusing 



lens, the 


distance between 


the 


focusing 


lens 


and 


the 


particles , 


the 


diameter of 


the 


motor window. 


and 


the 


height/position 


of the diode 


array 


[Ref. 1]. 


The 


range 


of 
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scattering angles for this investigation was .01-. 06 
radians, corresponding to an offset of 5 nun of the array and 
.02-. 07 radians corresponding to a 10 nun offset. Because of 
the very small detectors, the system is not capable of 
measuring very weak scattering intensities. 

The second system utilized was a commercially available 
particle size analyzer manufactured by MALVERN Instruments. 
The particular instrument used in this series of experiments 
was the MALVERN 2600 HSD. Figure 3.7 depicts the MALVERN 
2600 system, which utilizes a 2 mW He-Ne laser with a 
collimated beam diameter of 9 mm (Table 3.2 details laser 
specifications) . One of the principle advantages of the 
MALVERN system is its use of an annular photodiode array. 
The array is a monolithic detector which consists of 31 
rings. It has increasing thicknesses of the annular detec- 
tor elements with increasing circumference, resulting in an 
increase of detector area as radius increases. Ring arrays 
permit much weaker intensities to be measured than can be 
accomplished with linear arrays. The dynamic range of the 
instrument is 180:1 on any of the three range settings, with 
an advertized accuracy of plus or minus 4% of volume mean 
diameter [Ref. 19]. The detector may be used in conjunction 
with any one of three interchangeable Fourier transform 
lens, with focal lengths of 63, 100 and 300 mm [Ref. 14]. 
The respective size classes associated with each lens are 
1.2-118, 1.9-188 and 5.8-564 microns, with a sub-class 
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Figure 3.5 NFS Light Scattering Apparatus 





Figure 3.6 



Hewlett Packard HP 9836S Computer and 
HP 694 2A Multiprogrammer 
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Receiver Unit 




2600 Optical Unit - Components & Controls 




Figure 3.7 MALVERN 2600 Particle Sizer [Ref. 14] 
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TABLE 3.2 



LASER SPECIFICATIONS 



A. Helium-Neon Laser (Locally Designed System) 



B. 



1 . 


Manufacturer: 


Spectra-Physics 


2. 


Model : 


147 


3 . 


Type: 


He-ne Class 111B 


4 . 


Output Power: 


8 mWatt 


5. 


Beam Diameter: 


. 92 mm 


6. 


Beam Divergence: 


.87 mrad 


Helium-Neon Laser (MALVERN 


System) 


1 . 


Manufacturer : 


MALVERN Instruments 


2. 


Model : 


2600 


3. 


Type: 


He-Ne Class 111B 


4 . 


Output Power: 


2 mWatt 


5. 


Beam Diameter: 


9 mm 


6. 


Beam Divergence: 


Not Available 



which extends down to .5 microns [Ref. 14]. The 100 mm lens 
was used for this set of experiments. 

D. DATA ACQUISITION AND REDUCTION 

The data acquisition system associated with the locally 
designed laser diffraction system consisted of a Hewlett 
Packard 6942A multiprogrammer which was controlled by an HP 
9836S computer. The data acquisition program developed by 
Harris [Ref. 13] was the system driver in the acquisition 
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mode. In the course of a motor firing sequence, a 
background measurement was recorded prior to ignition. This 
mean intensity profile was then subtracted from the mean 
intensity profile measured with particles present. This 
procedure accounts for noise present due to ambient light 
and corrects for the characteristics of individual diodes. 
The raw data were displayed on the HP 9836S CRT so that any 
erroneous scans could be excluded prior to data reduction. 

A symmetric moving average type digital filter was then 
applied to the profile in order to smooth the curve [Ref. 
1] . This type filter does not have the phase lag that is 
inherent in analog filters. The data reduction algorithm is 
based on the work done by Buchelle [Ref. 10], and is fully 
described in Reference 20. The chief disadvantage of the 
method lies in the fact that there is no way to account for 
the actual distribution of the particulates. The Sauter 
mean diameter (SMD) determined by this algorithm will be 
that which is associated with the dominant mode (within the 
field of view of the array) of a multimodal polydispersed 
sample . 

The MALVERN system represents a vast improvement, in 
that not only is SMD determined, but the volume distribution 
is characterized as well. The MALVERN 2600 data acquisition 
system is centered about an AT&T 6300 PC. The system driver 
is the MALVERN Master-Sizer operating system. This version 
of the MALVERN software uses a command-based language that 
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allows the user to specify the pertinent parameters 
associated with the operation of the instrument. The 
Master-Sizer program allows the user to select from any one 
of four models. The model types available are Rosin- 
Rammler, log-normal, normal, and model independent. The 
model independent option of the Master-Sizer program is a 
process of constrained least squares fitting of theoretical 
scattering characteristics to observed data [Ref. 14]. It 
is usual to allow the best fit result to be obtained with no 
assumed form of size distribution and this allows the 
characterization of polydispersed multimodal distributions 
with high resolution [Ref. 14]. The capability of 
characterizing multimodal distributions is invaluable in 
work dealing with the combustion of metallized propellants, 



in that the 


distribution 


of particulates 


is 


most 


often 


multimodal . 


It 


is 


for 


this reason 


that 


the 


model 


independent 


mode 


was 


used exclusively 


throughout 


this 



investigation . 



24 



IV. RESULTS AND DISCUSSION 



A. INTRODUCTION 

This investigation was based on three major objectives: 

(1) refine previously developed experimental and 
diagnostic techniques by identifying specific causal 
factors of measurement error, 

(2) incorporate the use of a commercial particle size 
analyzer (MALVERN 2600) to validate and quantify the 
operating limits of the locally designed laser 
diffraction system and, 

(3) use the apparatus to determine the effects of 

operating pressure and position within the motor and 
exhaust nozzle on particle size distribution 

associated with the combustion of a reduced smoke 
propellant . 

Previous investigations at NPS had been based on the method 
of Buchele [Ref. 10] which yields Sauter mean diameter, but 
does not provide any information as to the actual 

distribution of particulate size. The introduction of the 
MALVERN system allowed for a characterization of volume 
distribution. The common parameter between the locally 

designed system and the MALVERN was Sauter mean diameter 
(D 32 )- This parameter provided the basis for comparison of 
the two diffraction systems. Both systems were used in the 
initial evaluation of the reduced smoke propellant. Once 
the limitations of the locally designed instrument were 
established, the MALVERN was used more extensively to gather 
data in both the motor and exhaust. The distinct advantage 
of being able to characterize multimodal particle size 
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distributions made the MALVERN the better suited instrument 
for the through motor measurements. In the exhaust, the 
annular detector of the MALVERN system allowed for the 
measurement of much weaker intensities as opposed to the 
locally designed system (which is also referred to as the 
NPS system) . 

B. SYSTEM CALIBRATION 

Among the methods available for calibration of laser 
diffraction instruments, the more commonly used include the 
calibration reticle and use of commercially prepared samples 
of known diameter. The latter approach was used in this 
investigation. The types of samples used included 
polystyrene spheres, glass beads and aluminum oxide powders. 
In each case, the sample material was suspended in a 
solution of distilled water. 

The locally designed system was calibrated using the 
polystyrene spheres, with diameters of 5.1, 9.6 and 2 0.0 
microns, respectively. Table 4.1 contains the results of 
these calibration tests. Figures 4.1 and 4.2 are 
representative samples of calibration results. 

The MALVERN 2600 was calibrated using three different 
materials, polystyrene spheres, glass beads and an aluminum 
oxide powder. Results are presented in Table 4.1. 
Calibrations from the two systems were in good agreement. 
System alignment was completed with the container of water 
in the optical path and a background measurement was taken 
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NPS IASER DIFFRACTION SYSTEM 
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CURVE FIT RESULTS 
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Figure 4.1 Calibration Profile for a Polvdispersed Sample (14.9-24.4 Microns) 
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Fiyure 4.2 Calibration Profile for a Polydispersed Sample (8.5-10.4 Microns) 



prior to introducing particles. A glass rod was used to 
stir the mixture to ensure that the particulates remained in 
suspension during the measurement process. 

A primary concern early in the investigation was to 
validate the ability of the MALVERN to accurately 

characterize a multimodal distribution when the system was 
operated in the Model Independent mode. In order to test 
this facet of the Master-Sizer program, a polydispersed 
bimodal sample consisting of 5.1 and 20.0 micron particles 
was prepared. The results are depicted as a volume 
distribution in Figure 4.3. The respective modes are 
centered at approximately 5.5 and 21.0 microns. 

When multiple scattering becomes significant, the mean 
particle diameters determined from single scattering theory 
are too small. Gulder [Ref. 281 has determined empirical 
corrections for the Malvern when obscuration is greater than 
0.5. This correction was utilized, when required, in the 
Dresent investigation. 



C. RESULTS FROM REDUCED SMOKE PROPELLANT EVALUATION 
1 . Test Conditions 

The series of experiments conducted covered a range 
of pressures from approximately 70 to 450 psig. The 
experiments were divided into three pressure classes, 150, 
250 and 400 psig. At each class of chamber pressure. 
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Figure 4.3 Volume Distribution, Bimodal Calibration Sample 




scattering measurements were made at the following 
locations : 

(a) near grain (end burning) 

(b) prior to the nozzle entrance 

(c) exhaust (P e = P C ritical) 

(d) exhaust (P e = P Q ) . 

A total of 32 tests were completed, 24 of which produced 
usable data. Three tests failed to produce data because the 
burst disk ruptured, three tests of the short motor were 
invalidated because the nitrogen purge system failed to keep 
the window surfaces clean, one test failed because of an 0- 
ring failure in the motor joint, and for one test, data were 
not obtained because the pressure threshold to trigger the 
acquisition of data was set too high. 

The nozzles were designed/sized for isentropic flow 
conditions using y = 1.2. The ratio of specific heats was 
determined using PEPCODE , an equilibrium, adiabatic 
combustion code. The value of gamma obtained in this manner 
was 1.21. Because nozzle efficiencies should be quite high 
(97-99%) , the isentropic assumption was made. Deviations in 
pressures can be attributed to the variations in ignition 
characteristics and propellant properties from run to run, 
and some small amount of nozzle throat erosion. 

The propellant composition was provided by the Naval 
Weapons Center at China Lake, California. The zirconium 
carbide in the propellant varied in size from 2-20 microns, 
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with a D 32 of 5.1 microns. The information is provided here 



as an introduction to the discussion of surface 
agglomeration on the burning grain. 

2 . Zirconium Carbide Behavior 

The size of agglomerates that form on the surface of 
burning metallized solid propellants are determined by many 
factors, including metallic powder type, propellant type, 
size of original metal particulates, concentration of metal, 
motor chamber pressure, and the burning rate of the 
propellant. The effect of chamber pressure on agglomerate 
size has been observed to be quite strong for aluminum. As 
chamber pressure and burning rate increase, the agglomerates 
become smaller [Refs. 21,22]. 

Zirconium carbide has been observed to burn 
vigorously on the surface of reduced smoke propellants 
containing AP and HTPB [Ref. 23]. If zirconium carbide 
(ZrC) burns stoichiometrically with oxygen it results in a 
temperature of approximately 4400°K. This resultant 
temperature is much greater than the melting temperature of 
Zr0 2 (2988°K) . Table 4.2 delineates the temperatures 
relevant to the discussion. 

A plausible behavior of the ZrC is that as the 
material is exposed at the grain burning surface, oxygen 
(0 2 ) from the ammonium perchlorate (AP) flame attacks the 
carbon to produce either carbon monoxide (CO) or carbon 
dioxide (C0 2 ). Oxygen from water vapor may also be 
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TABLE 4.2 



THERMOPHYSICAL PROPERTIES OF Zr, ZrC, AND Zr0 2 



Material 

Zr 

ZrC 

Zr0 2 



Melting 

Temperature f°K) 
2127 [Ref. 24] 
3523 [Ref. 26] 
2988 [Ref. 24] 



NOTE: Adiabatic flame temperature for 

determined to be 2916°K and the 
temperature for the reaction 

ZrC + 20 2 Zr0 2 + C0 2 

is approximately 4400°K. 



Boiling 

Temperature (°K) 
4700 [Ref. 24] 



4573 [Ref. 25] 

the propellant was 
adiabatic flame 



involved. The oxygen then reacts with the zirconium to 
produce a porous Zr0 2 surface layer [Ref. 24]. The oxide 
should be a liquid since the local temperature at the 
particle surface (as a result of the surface combustion of 
ZrC or Zr with 0 2 ) should be much greater than the melting 
temperature of Zr0 2 (2988°K) . The low burning rate 
associated with the relatively low chamber pressures allows 
time for the molten Zr0 2 on adjacent particles to stick 
together to form large agglomerates on the grain burning 
surface, prior to passing into the flow. At higher burning 
rates (pressures) the particles would be ejected from the 
surface more rapidly, resulting in smaller agglomerates. 
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3 . Near-Propel lant Surface Data 



The test results obtained near the grain surface are 
summarized in Table 4.3. D 32 values are plotted in Figure 
4.4. As can be seen directly from these results, there was 
a strong dependence on chamber pressure of the agglomerate 
size produced from the ZrC. At the lower pressures (159 and 
188 psig) , the distribution of particulates ejected from the 
grain surface were trimodal and bimodal, respectively. 
Figures 4.5 and 4.6 detail the characterization of the 
particulate distributions. The dominant modes were centered 
at approximately 30-37 microns and 65 microns. 

The Sauter mean diameters determined from the 
scattering data were 63.5 microns for the run at 159 psig 
and 72.8 microns for the run at 188 psig. This was 
considered to be good repeatability, run-to-run, allowing 
for variations in the combustion process from grain to 
grain. As discussed above, many of the ZrC particles, upon 
reaching the burning surface, apparently react to form Zr0 2 
and agglomerate before being ejected into the gas flow. 

Test number three of the near grain results (248 
psig) showed a reduction in Sauter mean diameter, with a 
measured D 32 of 32.2 microns. The distribution of 
particulates at this pressure was still essentially bimodal 
as depicted by Figure 4.7. However, the dominant modes of 
the distribution were now centered at approximately 12-13 
microns and 27 microns. There was also a small number of 
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Figure 4.5 Volume Distribution Near Grain Test No. 1 (159 psig) 
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Figure 4.6 Volume Distribution, Near Grain Test No. 
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Figure 4.7 Volume Distribution, Near Grain Test No. 3 (248 psig) 



particulates in the range of 3-4 microns, indicating that 
some of the particulates at this pressure (248 psig) were 
ejected directly into the gas stream. The 448 psig test 
resulted in a D 32 of 15.2 microns. Figure 4.8 reveals that 
the distribution of particulates continued to be bimodal. 
It is apparent from Figure 4.8 that a much larger portion of 
the original metal (ZrC) particles was ejected directly into 
the gas flow, as evidenced by the well-defined mode centered 
between 3-4 microns. The general observation here was that 
as motor chamber pressure was increased, the size of 
agglomerates on the grain surface decreased, similar to the 
behavior of aluminized propellants reported in References 21 
and 22. It was also observed that a greater number of the 
original ZrC particulates were apparently immediately 
ejected into the gas flow as they reached the burning 
surface as chamber pressure was increased. In all cases, 
the distributions of particulates near the burning surface 
were bimodal or trimodal. The quantities of particles in 
the smallest modes of the trimodal data were quite small. 

The most difficult obstacles to overcome in 
obtaining near grain light scattering measurements were: 
( 1 ) keeping the window surfaces clean, as the ignition 
process (due to its close proximity to the windows) made it 
difficult to effectively purge these surfaces, and ( 2 ) 
multiple scattering effects due to dense concentrations of 
particulate matter in the near grain region. The problem of 
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Figure 4.8 Volume Distribution Near Grain Test No. 4 (448 psig) 



multiple scattering has been addressed above. The empirical 
corrections presented by Gulder [Ref. 28 ] are for a Malvern 
system operated in the model independent mode. This mode was 
also used in the current investigation. When obscuration is 
approximately 0.99, the measured is approximately a factor 

of two too small. Gomi [Ref. 27 ] has also discussed the effects 
of multiple scattering. For monodispersed spheres Netzer and 
Powers [Ref. 1 ] have also presented data for correcting 
for high obscurations. 

The more difficult problem that had to be overcome 
in the near grain region was effectively purging the window 
surfaces. The only alternative here was to repeat the 
experiment until usable data were obtained. As noted 
earlier, three of the tests conducted in the near grain 
region were invalidated due to the fact that the window 
surfaces had become contaminated as a result of the ignition 
process . 

4 . Nozzle Entrance Data 

Table 4.4 presents the results from the measurements 
obtained prior to the nozzle entrance. These data allow one 
to examine the change in particle size across the nozzle 
when compared to the exhaust data which are contained in 
Tables 4.5 and 4.6. The D 32 data obtained using the MALVERN 
are also plotted in Figure 4.16. Figure 4.17 plots these 
same data together with the MALVERN data obtained near the 
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MOTOR CHAMBER SCATTERING MEASUREMENTS OBTAINED UPSTREAM 
FROM NOZZLE CONVERGING SECTION 
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es data obtained by Pruitt [Ref. 17]. 

es data obtained from NPS laser diffraction system; all other data 
obtained using the MALVERN 2600. 



burning surface. Figures 4.9, 4.10 and 4.11 characterize 
the volume distributions associated with the 150 psig class 
of experiments. Note in Figures 4.9 and 4.10 that the dis- 
tributions remained bimodal in nature, with the dominant 
modes being centered at approximately 35 and 68-70 microns 
in the case of Figure 4.9 and at approximately 27 and again 
at 68-70 microns in the case of Figure 4.10. In the case 
of Figure 4.11, note that the profile was trimodal, but 
shifted toward the smaller end, with dominant modes centered 
at approximately 3-4 microns, 20 microns and at 38-40 mi- 
crons. The Sauter mean diameter obtained from the scatter- 
ing data was 30.7 microns (Test No. 3). A post-test 
alignment check done with the MALVERN system in the align 
mode revealed that some significant window contamination had 
occurred, but physical examination of the windows showed 
very little evidence of contamination. The results of this 
one test, although questionable, are included for 
completeness. Comparison of the Sauter mean diameters and 
volume distribution profiles for near grain and nozzle 
entrance measurements completed for the 150 psig class 
reveal only slight differences. This occurred even though 
the residence time (t res ) was approximately 45 milliseconds. 

As discussed above the equilibrium adiabatic flame 
temperature was 2916°K. The actual temperature may have 
been somewhat less than this due primarily to heat loss 
through the uninsulated, stainless steel walls of the 
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Figure 4.9 Volume Distribution Nozzle Entrance Test No. 1 (119 psig) 
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Figure 4.10 Volume Distribution, Nozzle Entrance Test No. 




48 



Figure 4.11 Volume Distribution, Nozzle Entrance Test No. 3 (161 psig) 



combustor, and to the nitrogen gas used to purge the viewing 
windows. The observed particle behavior would be expected 
if the ZrC is completely consumed near (or on) the 
propellant surface, resulting in liquid Zr0 2 particles 
flowing in a low velocity (M << .1) environment. 

The NPS scattering data included in Table 4.4 will 
be discussed later in conjunction with the comparison of the 
MALVERN and NPS measurement systems. 

Figures 4.12 and 4.13 represent the volume distribu- 
tions associated with the runs completed for the 250 psig 
class. This represents data for tests 8 and 9 of Table 4.4. 
The Sauter mean diameter (SMD) determined from scattering 
measurements was 49.0 microns corresponding to test number 8 
and 48.0 microns associated with test number 9. This again 
represented good test-to-test repeatability. The dominant 
modes were centered at approximately the same locations for 
both tests. Comparison of these data with that obtained in 
the near grain surface region showed that there was approxi- 
mately a 16 micron particle growth over the length of the 
motor. Note that only one measurement was taken in the near 
grain region. If the point was valid, then particle 
collisions and/or wall collisions/ shedding may be suspected. 

Test numbers 10 and 11 represent data obtained for 
the 400 psig class runs. The SMD's determined were 32.0 and 
27.6 microns, respectively. Figures 4.14 and 4.15 are the 
volume distribution profiles for tests 10 and 11, 
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Figure 4.12 Volume Distribution, Nozzle Entrance Test No. 8 (252 psig) 
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Figure 4.13 Volume Distribution, Nozzle Entrance Test No. 9 (251 psig) 
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Figure 4.14 Volume Distribution, Nozzle Entrance Test No. 10 (370 psig) 
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Figure 4.15 Volume Distribution, Nozzle Entrance Test No. 11 (427 psig) 



respectively. Note that each was monomodal in nature. By 
comparing these profiles to those obtained near the grain 
surface, it was observed that the 3-4 micron diameter 
particles were eliminated as they passed along the length 
of the combustor. Again, if ZrC is assumed to react to ZrO ^ 
on or near the propellant surface, particle collisions could 
explain this observation. There was an appreciable increase 
in particle diameter from near the grain surface to the 
nozzle entrance section. Figure 4.16 is a graphical depiction 
of SMD plotted against chamber pressure. Note that the SMD 
decreases with increasing chamber pressure. Figure 4.17 
presents the data for both near-grain and nozzle-entrance 
measurements. From this figure, it is clear that particle 
growth/agglomeration occurred across the length of the 
combustor for pressures above 200 psi. 

5 . Exhaust Nozzle Data 

Table 4.5 is a listing of test results obtained in 
the exhaust for the condition of nozzle exit plane pressure 
equal to the critical pressure at the nozzle throat. The 
nozzles used in these tests had a 45° converging section and 
a short flat (.01"-. 015") in the throat section. All nozzles 
were truncated at the throat after the .01"-. 015" flat. 

Figures 4.18, 4.19, 4.20, 4.21 and 4.22 present the volume 
distributions associated with tests one through five, respect- 
viely. Figure 4.23 presents the results obtained on 
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Figure 4.18 Volume Distribution, Exhaust (P e = Pcritical^ Test No. (147 psig) 
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Figure 4.19 Volume Distribution, Exhaust (P e = Pcritical^ Test No. 
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Figure 4.20 Volume Distribution, Exhaust (P = P . . . -,) Test No. 3 (450 psig) 
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Figure 4.22 Volume Distribution, Exhaust (P = P . . • , ) Test No. 5 (371 psig) 

6 C IT 1 L 1 C d -L 



CURVE FIT RESULTS 
INTENSITY vs.THETR 




3DU110A QSdSiniJ 



63 



Figure 4.23 Exhaust Profile (P = P . . Test No. 6 (370 psig) 

L C -I. J. I — L C a -L 



test 6 using the NPS system. The distance downstream of the 
nozzle exit plane varied slightly as noted in Table 4.5 due 
to the geometry of the apparatus and set-up. All the 
distributions were monomodal with the exception of the test 
conducted at the lowest pressure (147 psig) . Here a second 
broadly distributed mode was observed, which covered a range 
of approximately 4.5-12.0 microns. Except for this point, 
the SMD varied only slightly as chamber pressure was increased. 

The Weber number is defined as 

w e = ^ag Qg ( u g “ u ag) 0 



where : 



d a g = the diameter of an agglomerate (m) 

pg = the density of the gas 

u a g = the velocity of the agglomerate 

Ug = the velocity of the gas 

0 = the surface tension. 

It is the ratio of inertial forces to surface tension 
forces. Gany [Ref. 22] found the Weber number to be a good 
correlation parameter in the analysis of the break-up of 
aluminum oxide agglomerates in the nozzle environment. Gany 
[Ref. 22] notes that for low Weber numbers the droplets are 
spherical, and tend to distort when the Weber number exceeds 
about four (for aluminum oxide). The distortion continues 
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to increase as Weber number increases until break-up occurs 
for Weber numbers in the range of 12-20 [Ref. 22]. 

Figure 4.24 is a plot of SMD versus motor chamber 
pressure. It is apparent from the curve that the 

agglomerates broke up into nearly the same size in the 
converging section of the nozzle, irrespective of chamber 
pressure (unless very low) or particle size entering the nozzle. 

The equilibrium adiabatic combustion temperature was 
approximately the same as the melting temperature for Zr0 2 . 
Thus, even though some thermal lag can be expected, the 
particles could be expected to rapidly solidify as the gas 
expansion in the nozzle occurs. This would prevent any 
further reduction in particle size. 

Table 4.6 contains the results obtained in the 
exhaust for the condition of nozzle exit pressure equal to 
ambient pressure. The nozzles all had converging half 
angles of 45° and diverging half angles of 15°. Nozzle 
throat diameter and exit diameter were determined using 
isentropic relations for a ratio of specific heats equal to 
1.2. Figures 4.25 to 4.32 present the volume distributions 
associated with tests 1 through 8, respectively. It is 

interesting to note that all the distributions were 
multimodal with the exception of test number 5, which was 
run at 215 psig. The primary factor which differentiated 
test five from the other tests was the distance behind the 
nozzle exit plane at which the measurements were obtained. 
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EXHAUST SCATTERING MEASUREMENTS (P 
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NOTE: All data were obtained using the MALVERN 2600. 
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Figure 4.25 Volume Distribution, Exhaust (P = P ) Test No. 
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Figure 4.26 Volume Distribution, Exhaust (P = P ) Test No. 
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Figure 4.27 Volume Distribution, Exhaust (P = P ) Test No. 3 (135 psig) 
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Figure 4.28 Volume Distribution, Exhaust (P^ = P ) Test No. 4 (79 psig) 
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Figure 4.29 Volume Distribution, Exhaust (P = P ) Test No. 5 (215 psig) 
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Figure 4.30 Volume Distribution, Exhaust (P = P ) Test No. 6 (436 psig) 
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Figure 4.31 Volume Distribution, Exhaust (P e = P Q ) Test No. 7 (352 psig) 




75 



Figure 4.32 Volume Distribution, Exhaust (P = P ) Test No. 8 (300 psig) 



The measurements were all made 2.7 inches downstream of the 
nozzle exhaust plane except for test five, where the 
distance was 0.6 inches. The tests run with the 2.7 inch 
distance were necessary because laser diffraction data were 
being obtained simultaneously through the motor, and the 
physical dimensions of the apparatus limited the minimum 
distance to 2.7 inches. Test 5 results indicated that the 
majority of the particles remained monomodal and essentially 
unchanged in mean diameter (19 microns) from the throat to 
the exit plane. The measurements made further aft had one 
mode centered at the same 19 micron diameter. However, the 
number of these particles were less, and significant numbers 
of smaller particles were observed in a bi- or tri-modal 
distribution. The reduced particle size could be due to particle 
breakup. However, it may have resulted from the reduced 
number of particles in the scattering volume, which was locat- 
ed through the center of the spreading exhaust jet. 

The presence of strong thermal gradients along the 
boundary of the exhaust jet initially caused some problem in 
obtaining exhaust measurements. During initial runs, 

scattering data indicated the presence of abnormally large 
particulates in the exhaust, on the order of 127 microns. 
Knowing that particle size entering the nozzle was 
substantially less than 127 microns, typically around 70 
microns for 150 psig runs, a search for causal factors was 
pursued. Conversation with Hughes and Taylor at MALVERN 
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Instruments confirmed the idea that thermal gradients or 
density gradients were indeed biasing the scattered light 
data, i.e., the presence of strong thermal gradients along 
the boundary of the hot gas jet, or density gradients 
associated with a shock. The solution to this difficulty 
was a rather simple procedure. When unstable refractive 
index of the suspending medium is present as a result of 
strong thermal gradients, a steering effect occurs that 
steers light to low angles, hence biasing particle size to 
the large end of the scale. The procedure employed when 
using the Master Sizer program was to inhibit data from the 
affected inner rings, using the "KILLDATA" command. Once 
the user had examined the raw data and observed which rings 
had been affected, then those rings were inhibited by the 
software command "KILLDATA." Care must be exercised in this 
procedure to ensure that no scattering data will be lost as 
a result of inhibiting data from the affected inner rings. 
In the latter case there is no recourse, scattering data at 
that point becomes no longer separable from the steered 
light. 

Figure 4.33 is a graphical depiction of particle 
size (SMD) versus motor chamber pressure. There was a 
slight reduction in particle size with increasing pressure. 
Figure 4.34 plots both exhaust conditions. The reduced size 
in the downstream, spreading plume is apparent. 
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Figure 4.34 Comparative Results, Exhaust (P = p cr i t i ca i) and 
Exhaust (P = P ) Data 



Exhaust measurements were made for both P~= p and p„ = 

c O “ 

^critical to determine if any change in particle size was 
incurred as a result of the supersonic expansion of the jet. 
Particle breakup downstream of the nozzle exhaust was observed. 
The limited data at approximately 200 psi indicated that little 
breakup occurred within the nozzle, downstream of the throat. 

In aluminized propellants, one might reasonably expect that 
vapor condensation occurring during the expansion process might 
well affect the size distribution of exhaust particles. How- 
ever, in the case of zirconium carbide, the temperature in the 
chamber is not sufficiently high for appreciable amounts of 
either Zr, ZrC or ZrC >2 to be present in a vapor phase. 

6 . SEM Analysis 

As a means of independently verifying the scattering 
results, residue was collected from the face of the nozzle 
converging section. This residue was then placed in a solution 
of lab grade acetone and put into an ultrasonic cleaner, mixed, 
and allowed to settle. After repeated cleanings, the acetone 
was evaporated off in an oven under reduced pressure. Using 
this powder, seven samples were prepared for examination with 
the scanning electron microscope. Three of the samples were 
prepared using a non-conducting epoxy base which had to have 
a gold flashing applied in order to provide for a conducting 
surface. The remaining samples were mounted using a conduct- 
ing colloidal carbon. Figures 4.35 and 4.36 are SEM photo- 
graphs taken of the samples. Particles were counted from each 
sample in an attempt to get a representative size distribution. Seme 218 
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Figure 4.35 SEM Photograph of Residue Collected 
from Exhaust Nozzle 
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Figure 4.36 SEM Photograph of Residue Collected 
from Exhaust Nozzle 
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particles were counted with sizes ranging from 1-20 microns. 
The smaller particles (3 microns and less) were difficult to 
see, but did exist in the sample. No particle larger than 
20 microns was observed, but this by no means implies that 
no larger particles were present. Collection of residues 
and/or exhaust products is a crude method. When collected 
in small quantities from wall surfaces it cannot yield a 
valid statistical sample. However, it does provide some 
reassuring insight as to what was physically present in the 
sample. 

The particles were grouped strongly in 2 bands, one 
band from 5-12 microns and the other from 14-17 microns. 
The sample was taken from the 427 psig run where the mean 
particle size (SMD) prior to the nozzle entrance was 27.6 
microns and in the nozzle throat region it was 19 
microns. Thus, the SEM data were in reasonably good 
agreement with the optical data. In Figures 4.35 and 4.3 6 
the larger irregular-shaped particles have reasonable 
sources, one being residue from the binder/ inhibitor , and 
the second being residue from the lab grade acetone used to 
clean the samples. 

D. COMPARATIVE EVALUATION OF THE MALVERN 2600 AND THE 

LOCALLY DESIGNED LASER DIFFRACTION SYSTEM 

From the outset of this investigation, a comparative 
analysis of the two diffraction instruments was a primary 
goal. Some of the limitations of the NPS system were 
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obvious, based on a review of the system hardware and 
software . 

The differences in software and in data reduction 
algorithms were significant. The NPS algorithm is based on 
a ratioing technique developed by Buchele [Ref. 10]. This 
method yields the size parameter D 32 by the process of 
fitting a theoretical curve for monomodal distributions to a 
measured intensity profile. The technique is discussed in 
detail in Chapters II and III. The MALVERN Master-Sizer 
operating system, by comparison, inverts the scattered data 
to obtain a size distribution. The NPS system could be 
upgraded in this area by adopting one of the quadrature 
methods described by Koo [Ref. 16] . 

The differences in hardware are highlighted by the 
design of the detector system photodiode arrays. The 
MALVERN systems annular array design allows for the 
detection of much weaker intensities. The field of view of 
the NPS system was also limited by the array dimensions. 

A comparison was made using the NPS laser system to 
measure particle size at the nozzle entrance. Figure 4.37 
is a plot of comparable through-motor measurements made by 
the MALVERN system and the NPS system in the motor section 
just prior to the nozzle entrance. As can be seen from the 
plot, the NPS system underestimated the mean particle sizes 
at low pressures where the distributions were multimodal. 
As discussed above, the NPS data reduction algorithm assumes 



83 




(snohdiw) aaiawvia Nvaw nainvs 



O 

O 

lO 



0 ) 

x: 

4-) 



o 






o 


u 


5m 




H 


0 




u} 


UM 




Cu 






' — ' 


03 






-p 




w 


03 




(X 


Q 




D 






CO 






CO 


os 




W 


w 


o 


PS 


> 


o 


cu 


pi 


m 


ps 






w 






CQ 






£ 


c 






03 










u 


CO 






Oh 




OS 




o 


O 




o 


Eh 


«* 


CM 


o 


+J 




£ 


i—l 



D 



in 


c 


a 


c 


os 


c 




o: 


0 


u 


> 


4i 


•H 


C 


+J 


C4 


03 




U 


cu: 


03 


«— i- 


a 


N, 


£ 


N 


0 


o: 


u 


2 



m 

O <D 

U 

D 

•H 

Cm 



84 



a monomodal distribution. The larger particles in the bimodal 
distribution (60-70 microns) , scattered light at angles below 
the range of the NPS detector. To further verify this effect, 
the MALVERN data for the run at 161 psig was reduced, after 
inhibiting data that were outside the field of view of the NPS 
system. Figure 4.38 depicts the resulting distribution. What 
occurred was that an artificial mode appeared at 24.7 microns. 
Figure 4.39 is a measurement made at a similar pressure (175 
psig) with the NPS system. A SMD of 24 microns was obtained. 
Thus, when distributions are multimodal, accurate values of 

are not obtained with present data reduction algorithm 
used with the NPS apparatus. Figures 4.40, 4.41 and 4.42 are 
included for completeness. These figures represent the data 
used to construct Figure 4.37. At higher chamber pressures 
(400 psig class) , where the distributions became more mono- 
modal and the dominant scattered light was within the field 
of view of the NPS system, better agreement of data existed 
(Figure 4.37) . 

The conclusion is that the NPS system functions well when 
the distribution of particulates is monomodal and the linear 
diode array is placed such that the majority of the scattered 
light is within the field of view of the detector. 

E. SOURCES OF MEASUREMENT ERROR (UNCERTAINTY) 

In the course of this investigation a number of 
measurement errors were encountered from a variety of 
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Figure 4.38 Example of the Effects of a Reduced Field of View 
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Figure 4.39 Motor Profile, Nozzle Entrance Test No. 
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Figure 4.40 Motor Profile, Nozzle Entrance Test Mo. 5 (125 psig) 
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Figure 4.41 Motor Profile, Nozzle Entrance Test No. 6 (135 psig) 
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Figure 4.42 Motor Profile, Nozzle Entrance Test No. 7 (79 psig) 



sources . 



Measurements were made in both the motor chamber 



and exhaust environments, and each environment presented a 
slightly different set of difficulties. 

Measurements made in the motor environment were, subject 
to the following sources of uncertainty: 

(a) scattering by the aperture (window frame) 

(b) limited field of view 

(c) multiple scattering effects due to dense 

concentrations 

(d) surface reflections from the windows in the motor 

(e) diameter of aperture limited scattering angle to a 
maximum of 6° 

(f) cleanliness of window surfaces 

(g) extraneous light 

(h) quenching effects due to chamber walls, and the 
introduction of nitrogen used as a purge for the 
window surfaces 

(i) vignetting of scattered light, i.e., sample not 
completely within the range of the focal length of 
the receiving lens. In the event vignetting occurs, 
particle size is biased to the larger sizes. 

Scattering by the aperture is discussed in some detail 
by Buchele [Ref. 10]. Although some scattering from this 
source does occur, if the particle concentration is 
sufficiently high (as in the present case) , particle 
scattering will dominate. 

The limited field of view of the NPS system did 
significantly impact the scattering data obtained. In the 
case of the present investigation, the physical dimensions 
of the photodiode array did not allow for scattered light to 
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be detected at angles less than .01 radians. Also in the 
current configuration of the NPS scattering device, the largest 
particle size which could reasonably be detected was 40 microns 
The problem may be overcome in large part by increasing the 
dimension of the array. A 4096 diode RETICON array has been 
purchased, and should be installed in the near future. This 
will provide a much enhanced field of view for the NPS instru- 
ment . 

Multiple scattering effects tend to occur when the separat- 
ion distance between particles becomes small. Correction 
methods have appeared in the literature [Refs. 27,28]. These 
effects were present in the motor measurements and empirical 
corrections of the data were required. A reduced combustion 
chamber width is required to reduce this problem. 

The problems associated with windowed surfaces are three- 
fold in nature. The first is that some reflections are present 
The MALVERN system in the align mode allows sufficient adjust- 
ment of the detector in order to eliminate the effects due to 
reflections. The second problem was a result of the physical 
diameter of the window which could pass light. This restricted 
the maximum scattering angle to 6°. If very small particles 
had been present, this would have had a biasing effect on the 
received scattering data. The apparent size would have been 
larger than the actual size. The third problem associated with 
the windcwed surface was to provide sufficient purge to keep these 
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surfaces clean. A trial and error method of adjusting purge 
flow rates was employed. If the windows had become 
contaminated in a previous run, the flow rate was increased 
until the surfaces remained clean and usable data were 
obtained. 

The introduction of the purge nitrogen, along with the 
walls of the motor chamber, had some undetermined quenching 
effect on the gas and particle temperatures. Although the 
percent of gas flow comprised by the nitrogen purge was 
estimated to be between 3 and 9 percent for the given range 
of pressures, it was not possible to determine what 
influence the quenching effect had on particle size, or the 
extent to which it actually occurred. 

Extraneous light can also be a factor. Buchele [Ref. 
10] details these effects. The ratioing method adopted in 
the data reduction algorithm for the NPS system has the 
desirable feature of eliminating this effect. Reference 10 
points out the fact that extraneous light is most likely to 
be a problem at small angles. 

Vignetting effects occur when scattered light is not 
captured by the receiving lens. This tends to bias particle 
size toward the larger size bands. The way to overcome this 
difficulty is to ensure that the field of particulates is 
well within the focal length of the receiving lens, so that 
all of the light scattered by the sample will be captured by 
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the (Fourier transform) lens and focused onto the detector 
array. 

The measurement errors that were found to be associated 
with obtaining scattering data in the exhaust are listed 
below: 

(a) unstable refractive index of the suspending medium 
which result from thermal or density gradients 

(b) turbulence along the jet boundary 

(c) shocks associated with a jet which is not ideally 
expanded 

(d) the location at which the measurements are taken 
aft of the nozzle can affect the measured 
distribution 

(e) vignetting, similar to that described previously 

(f) extraneous light (also discussed previously in 
conjunction with motor measurement errors) . 

Thermal gradients play a significant role in measuring 
particle size in hot gas jets. Close to the exit plane of 
the exhaust, the jet was slightly smaller than the MALVERN 
laser beam diameter. As discussed above, scattering data 
were influenced by the thermal gradients along the jet 
boundary. Also as discussed earlier, the "KILLDATA" command 
can be used to inhibit data from the inner rings where the 
steering effect is most prominent. Buchele [Ref. 10] points 
out that turbulence compounds these effects. It was 
impossible to separate the effects of turbulence and thermal 
gradients on the beam steering. In either case, so long as 
the user is reasonably assured that no scattering data are 
present as a result of large particles in the exhaust flow, 
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the "KILLDATA" command may be used to inhibit data from the 
light steered to low angles. 

The presence of shocks in the exhaust jet which is not 
ideally expanded, is to be expected. The density gradients 
in the region of a shock can also adversely affect scattered 
data. If shocks are identified as a problem, a schlieren 
may be used to establish shock location so that on 
subsequent runs the apparatus can be positioned upstream of 
the shock. 

The location at which diffraction measurements are made 
in a supersonic exhaust can influence the resulting 
distribution of particulates which are observed. The 
further back in the jet, the less concentrated the particle 
field becomes. This may have biased the measurements to 
the smaller particles which were present. 

The effects of extraneous light and vignetting have been 
previously discussed and will not be discussed further, 
except to say that the measurement and subtracting of 
background light plays a role in reducing the effect of 
extraneous light. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The results of this investigation have demonstrated that 
laser diffraction methods can be an effective means for 
determining particle sizes in solid-propellant rocket 
motors. Measurements were successfully made in both the 
motor chamber and exhaust environments. Each environment 
presented a different set of difficulties. 

The MALVERN and the locally designed diffraction system 
produced data which correlated quite well when the particle 
distributions were monomodal, and the dominant scattered 
light was within the field of view of the NPS linear diode 
array. When the distributions were multimodal, the NPS 
system did not accurately measure the Sauter mean diameter. 

This problem was the result of both a limited field of 
view and the data reduction algorithm. The latter was based 
upon monomodal distributions. The system could be 
significantly improved by using larger arrays and by using 
data reduction techniques which are based on quadrature 
methods [Ref. 16]. 

The NPS system will always be limited to higher number 
density flows since the linear array only intercepts a very 
small fraction of the light scattered at a given angle. 

The MALVERN 2600 has proven to be well suited for use in 
solid propellant combustion experiments. The beam diameter, 
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however, is larger than optimum for the current application. 

The acquisition of a new beam expander to reduce the diameter 
of the collimated beam would enhance both in-situ motor meas- 
urements and exhaust measurements. 

A transition from the axisymmetric motor to a two-dimensional 
motor would improve through-motor measurements in the following 
ways : 

(a) advantage of reduced optical depth 

(b) reduction or elimination of multiple scattering effects 
that may be present in the motor chamber environment 
due to dense particle concentrations 

(c) increased range of scattering angles 

(d) allows operation to higher pressures. 

The motor data suggested that ZrC reacts to form ZrC^ on 
or near the propellant surface. Surface agglomerations occur- 
red and were observed to decrease in size as chamber pressure 
increased . 

At the low Mach numbers in the present motor the ZrC >2 
particles appeared to remain practically unchanged as they 
passed down the length of the motor port at very low pressures. 

At higher pressures (above 250 psig) , particle growth/agglom- 
eration occurred. 

At low pressures the ZrC^ particles were multimodal. At 
higher pressures ( greater than 350 psi) the distributions 
approached monomodal. 
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ZrC >2 particles broke up in the nozzle converging section 
to a monomodal distribution with a between 15 and 19 microns, 

independent of the particle size entering the nozzle. Limited 
data indicated that no appreciable change in particle size 
occurred in the supersonic nozzle flow. Aft of the nozzle, 
in the expanding plume, the particle size was significantly 
smaller (2-8 microns). This could be due to particle breakup, 
or to measurement location biasing. 

SEM analysis of particles collected form the converging 
nozzle wall showed particles with sizes in good agreement 
with the optically obtained data. 
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